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We studied self-diffusion of small 2D Ni islands (consisting of up to 10 atoms) on Ni (111) surface 
using a self-learning kinetic Monte Carlo (SLKMC-II) method with an improved pattern-recognition 
scheme that allows inclusion of both fee and hep sites in the simulations. In an SLKMC simulation, 
a database holds information about the local neighborhood of an atom and associated processes that 
is accumulated on-the-fly as the simulation proceeds. In this study, these diffusion processes were 
identified using the drag method, and their activation barriers calculated using a semi-empirical 
interaction potential based on the embedded-atom method. Although a variety of concerted, multi- 
atom and single-atom processes were automatically revealed in our simulations, we found that these 
small islands diffuse primarily via concerted diffusion processes. We report diffusion coefficients for 
each island size at various tepmratures, the effective energy barrier for islands of each size and the 
processes most responsible for diffusion of islands of various sizes, including concerted and multi- 
atom processes that are not accessible under SLKMC-I or in short time-scale MD simulations. 

PACS numbers: 68. 35. Fx, 68.43.Jk,81.15.Aa,68.37.-d 



I. INTRODUCTION 

Surface diffusion is of interest not only because it is so 
different from diffusion in bulk solids 1 but because dif- 
fusion of adatoms on metal surfaces, individually or as 
a group via multi-atom or concerted diffusion processes 
plays an essential role in a wide variety of such surface 
phenomena as heterogenous catalysis, epitaxial crystal 
growth, surface reconstruction, phase transitions, segre- 
gation, and sintering. 2 A precise knowledge of diffusion 
mechanisms is essential for understanding and control of 
these phenomena. 3 Adatoms can diffuse on a substrate 
in a variety of ways, and competition between various 
types of diffusion processes (due to the differences in their 
rates) determines the shapes of the islands formed and 
(on macroscopic times scales), the morphological evo- 
lution of thin films. Hence a great deal of effort has 
been devoted to investigation of self-diffusion of adatom 
islands on metal surfaces, initially using field ion mi- 
croscopy (FIM), 4-9 and more recently scanning tunnel- 
ing microscopy (STM). 10-18 Because of inherent differ- 
ences in the microscopic processes responsible for island 
diffusion on different metal surfaces, this is still an on- 
going research problem. Both experimental and theoret- 
ical studies for various systems have succeeded in finding 
the activation barriers and prefactors for a single-adatom 
diffusion processes 19-31 . Ref 1 provides a good survey of 
those efforts. To the best of our knowledge, however, 
there has so far been no systematic experimental or the- 
oretical effort to identify the diffusion mechanisms re- 
sponsible for diffusion of small 2D Ni islands on Ni(lll) 
and to calculate their activation barriers. In this article 
we report our results of doing so for such islands, ranging 
in size from 1 to 10 atoms. 

Arrangement of atoms in the substrate of an fcc(lll) 
surface results in two types of three-fold hollow sites for 
an adatom: the regular fee site (with no atom beneath 
it in the second layer), and an hep site (with an atom 



beneath it in the second layer). Occupancy of adatoms 
at fee sites maintains the crystal stacking order (ABC 
stacking) of fee structure, while occupancy of hep sites 
leads to a stacking fault. Depending on its relative occu- 
pation energy, which is material dependent, an adatom 
can occupy one or the other of these sites. Which site is 
preferred on the fcc(lll) surface affects the way diffusion 
and hence growth proceeds. It is therefore important to 
understand whether the diffusion proceeds via movement 
of atoms from fcc-to-fcc or hep-to-hep or fec-to-hep hep- 
to-fee sites. It has been observed experimetally that for 
smaller clusters mixed occupancy 32 of fee & hep sites is 
possible. 

A host of studies has been devoted to problems of self- 
diffusion and diffusion mechanisms on metal fcc(lll) sur- 
faces, almost exclusively, however, with either a precon- 
ceived set of processes or merely approximate activation 
barriers. It is nevertheless crucial to discover the full 
range of processes at work and to accurately establish 
the activation barrier of each. It is also well known that 
the fcc(lll) surface, being atomically flat, has the least 
corrugated potential energy surface of any fee surface, re- 
sulting in low diffusion barriers even for clusters to diffuse 
as a whole. Consequently, studies of diffusion processes 
on fcc(lll) surfaces is a challenging problem for both ex- 
periment and simulation even to this day. For a monomer 
and smaller islands like dimer, trimer and up to certain 
extent, tetramer, all possible diffusion processes may be 
guessed. But as islands further increase in size, it be- 
comes more difficult to enumerate all possible diffusion 
processes a priori. An alternative is to resort to molec- 
ular dynamics (MD) simulation. But because diffusion 
processes are rare events, an MD simulation cannot cap- 
ture every microscopic process possible, as most of the 
computational time is spent in simulating atomic vibra- 
tion of atoms. Instead, to do a systematic study of small 
Ni island diffusion on Ni (111) surface we resorted to 
an on-lattice self-learning kinetic Monte Carlo (SLKMC- 
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II) method, which enables us to study longer time-scales 
than are feasible with MD yet to find all the relevant 
atomic processes and their activation barriers on-the-fly, 
as KMC methods limited to a priori set of processes can- 
not do. Moreover, whereas previous studies have used an 
on-lattice SLKMC method, 33-36 in which adatoms were 
restricted to fee occupancy, in the present study both 
fee and hep occupancies are allowed, and are detectible 
by our recently developed improved pattern-recognition 
scheme 37 . 

The remainder of the paper is organized as follows. In 
Section II we discuss the details of our SLKMC-II simu- 
lations, with particular attention to the way we find dif- 
fusion processes and calculate their activation barriers. 
In Section III we present details of concerted, important 
multi-atom and single-atom diffusion processes responsi- 
ble for the diffusion of Ni islands as a function of island 
size. In Section IV we present a quantitative analysis of 
diffusion coefficients at various temperatures and of ef- 
fective energy barriers as a function of island size. In 
Section V we present our conclusions. 



II. SIMULATION DETAILS 

To study Ni island diffusion on fee Ni(lll) surface, 
we carried out SLKMC simulations using the pattern- 
recognition scheme we developed recently 37 that includes 
both fee and hep sites in the identification of an atom's 
neighborhood. Various types of diffusion processes are 
possible, and their activation barrier depends on the 
atom's local neighborhood. Whenever a new neighbor- 
hood around an atom is identified, a saddle-point search 
is carried out to find all the possible atomic processes 
and calculate their activation barriers - provided that it 
has at least one similar empty site in the second ring 3 ' , 
since when an atom occupies an fee (or alternatively an 
hep) site, the nearest neighbor (NN) hep (or, correspond- 
ingly, fee) sites cannot be occupied. In our simulations 
we used a system size of 16x16x5 with the bottom 2 lay- 
ers fixed, and carried out saddle-point searches using the 
drag method. In this method a central or active atom 
is dragged in small steps towards a probable final posi- 
tion. If the central atom is on an fee (hep) site, then it is 
dragged towards a NN vacant fee (hep) site in the second 
ring. Since atoms are allowed to occupy either hep or fee 
sites, an atom being dragged from an fee (hep) site to a 
neighboring similar site is allowed to relax to an interme- 
diate hep (fee) site in between the two fee (hep) sites. In 
other words, processes are possible in which atoms in an 
island may occupy fee, hep or both fee & hep sites simul- 
taneously. In the drag method, the atom being dragged 
is always constrained in the direction of the reaction co- 
ordinate but allowed to relax along its other degrees of 
freedom (those perpendicular to the reaction coordinate) , 
while all the other atoms in the system are allowed to re- 
lax in all degrees of freedom. Once the transition state 
is found, the entire system is completely relaxed to find 



the final state of the process. The activation barrier of 
the process is the difference between the energies of the 
transition and initial states. We have verified the ac- 
tivation barriers of some of the key processes found by 
the drag method using the (more accurate but computa- 
tionaly expensive) nudged-elastic band (NEB) method 38 , 
and found no significant difference. For inter-atomic in- 
teractions, we used an interaction potential based on the 
embedded-atom method (EAM) as developed by Foiles 
et al 39 . In all our SLKMC simulations we used the same 
pre-exponential factor of 10 12 s _1 , which has been demon- 
strated to be a good assumption for such systems as the 
one under examination here. 40 ' 41 

For the small islands under study here (1-10 atoms), 
we found that when an atom is dragged rest of the atoms 
in the island usually follow. For very small islands (1-4 
atoms), all of the processes identified by the drag method 
were concerted-diffusion processes. As island size in- 
creases we found single-atom and multi-atom processes 
as well. For islands of size 5-6, even single-atom detach- 
ment processes are identified and stored in the database 
(even though they are not allowed in our simulations). 
To account for all types of processes associated with both 
compact and non-compact shapes - especially concerted 
processes and multi-atom processes - we used 10 rings to 
identify the neighborhood around an active atom in our 
SLKMC simulations. Using 10 rings corresponds to in- 
cluding fifth nearest-neighbor interactions. To make sure 
we identified all the single-atom processes, we also car- 
ried out saddle-point searches with all of the atoms fixed 
except the atom being dragged. Although there is no in- 
fallible method for discovering all possible processes, we 
did exhaust the search for possible processes identifiable 
using the drag method. 

In order to save computational time, we first carried 
out SLKMC simulations at 700K for each island size, 
and used the database thus generated to carry out our 
simulations for the same size at lower tempratures (300, 
400, 500 and 600K ). The rationale for this approach is 
that an island goes through many more shapes at higher 
tempratures: when a simulation is carried out at a lower 
temprature starting out with a database generated at a 
higher temprature, it only rarely finds an unknown con- 
figuration. It is not possible, however, to economize on 
computational time by using, for the smaller islands un- 
der study here, a database generated for (say) the larger 
among them, because the types of processes possible 
(along with their respective barriers) are dependent on 
an island's particular size. 



III. RESULTS 

As mentioned above, all of the processes for a given 
island are identified and their activation barriers calcu- 
lated, and stored in a database on-the-fly. We discuss in 
this section, however, only key processes of the various 
general types (concerted, multi-atom and single-atom). 
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FIG. 1. (a) fee and hep sites on an fcc(lll) surface, with 
corresponding directions for concerted diffusion processes; (b) 
A-type and B-type step edges (here, for an all-hep island) for 
the same surface. 



Fig. 1(a) is a sketch of the fcc(lll) surface with its 
adsorption sites marked as fee and hep. Determining 
whether an adatom is on an fee or on an hep site on this 
surface requires knowledge of at least 2 substrate layers 
below the adatom layer. In all our figures we show only 
the adatom layer and the layer below (the top substrate 
layer) with the convention that the center of an upward- 
pointing triangle (along the y-axis) formed by the (top 
layer) substrate atoms is an fee site, while the center of 
a downward-pointing triangle pointing triangle is an hep 
site. An island on an fee (111) surface can be on fee 
sites or on hep sites or a combination of both sites (some 
atoms of the island sitting on fee sites and the rest on 
hep sites) . Depending on the type of material either the 
fee or the hep site will be energetically favorable. As we 
shall see for each island size under study here, the fee 
site for Ni(lll) is always at least slightly more favorable 
than the hep site. 

A compact adatom island on an fcc(lll) surface can 
move in the three directions shown in Fig. 1(a). Note that 
the numbering scheme for the directions open to an atom 
on an fee site is inverse to that for those open to an atom 
on an hep site (see Fig. 1(a)). We follow the enumera- 
tion convention for directions distinguished in Fig. 1(a) 
throughout the article in tabulating activation barriers 
for concerted processes for islands of various sizes and 
shapes. Concerted processes involve all atoms moving 
together from all-fee sites to all-hep sites or vice-versa. 
In a concerted diffusion process a cluster can either trans- 
late in one of the three directions shown in Fig. 1 (con- 
certed translation) or rotate around an axis (around the 
center of mass), either clockwise or anti-clockwise (con- 
certed rotation). Since concerted rotational processes do 



not produce any displacement in the center of mass of 
an island, they do not contribute to island diffusion. De- 
pending on the size of the island and its shape, activation 
barriers for the processes in these three directions can be 
different. 

Activation barriers for single-atom processes, however, 
depend on the type of step-edge along which atom dif- 
fuses. Fig. 1(b) shows, using the example of a 6-atom hep 
island, how an A-type step-edge — a (100) micro-step dif- 
fers from a B-type step-edge — a (111) micro-step. We 
discuss important single-atom diffusion processes system- 
atically and in detail in Sub-section IIIK. 

As island size increases not only does the frequency 
of single-atom processes increase but the frequency of 
multi-atom processes does so as well. All multi-atom 
mechanisms involve shearing. A special case is reptation 
mechanism 42,43 , a two-step shearing process that moves 
the cluster from all-fee to all-hep sites or the reverse: 
first, part of the island moves from fee to hep sites; then 
the rest of the island moves from fee to hep. Hence at 
the intermediate stage, the island has mixed fec-hep oc- 
cupancy. In case of Ni-island diffusion, reptation pro- 
cesses occur only when the shape of the island becomes 
non-compact. We will discuss reptation in detail when 
we take up islands of size 8-10. 



A. Monomer 

As mentioned earlier, much work has been done to de- 
termine activation barrier for Ni monomer diffusion on 
Ni (111) surface 1 . A monomer on fee (111) surface can 
adsorb either on an fee or an hep site. We find that ad- 
sorption of an adatom on an fee site is slightly favored 
over than on an hep site by 0.002 eV - in good agree- 
ment with the value reported in Ref 22. Diffusion of a 
monomer occurs through hopping between fee sites via 
an intermediate hep site. We find the activation barrier 
for a monomer's hopping from an fee site to a neighboring 
hep site to be 0.059 eV while that for the reverse process 
is 0.057 eV. The effective energy barrier for monomer is 
found to be 0.057 eV, which is consistent with the result 
reported by Liu et al. 18 of 0.056 eV. 



B. Dimer 

On any fcc(lll) surface a dimer (of the same species) 
can have three possible arrangements: both atoms on fee 
sites (an FF-dimer, Fig. 2(b)), both on hep sites (an HH- 
dimer, Fig. 2(a)) or one atom on an fee and the other on 
an hep site (an FH-dimer, Fig. 2(c)). We find that the 
FF-dimer is energetically more favorable than the HH- 
dimer by 0.005eV and the FH-dimer the least favorable 
by 0.011 eV. We find that both FF and HH dimers diffuse 
via concerted as well as single-atom processes, whereas 
the FH-dimer diffuses via single-atom processes only. 
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FIG. 2. Possible configurations for a dimer, with activation 
barriers for concerted diffusion processes, (a) FF dimer (both 
atoms on fee sites); (b) HH dimer (both atoms on hep sites); 
(c) FH dimer (one atom on an fee and the other on an hep 
site); (d) FF dimer in concerted clockwise rotation and (e) 
HH dimer in concerted counter-clockwise rotation. 

TABLE I. Activation barriers (in eV) of concerted processes 
for dimer diffusion. 



Direction 


fee 


hep 


1 


0.071 


0.066 


2 


0.148 


0.143 


3 


0.148 


0.143 



In concerted diffusion processes, both atoms in a FF 
(HH) dimer move from fee (hep) to the nearest hep (fee) 
sites as shown in Fig. 2 (a)&(b) and Fig. 2 (d)&(e), 
thereby converting an FF (HH) dimer into an HH (FF) 
dimer. In the case of an FF (HH) dimer, the activa- 
tion barrier for concerted translational (Fig. 2 (a)&(b)) 
is 0.148 eV (0.143 eV) while that for concerted rotation 
(Fig. 2 (d)&(e)) is 0.038 eV. In concerted dimer rota- 
tion, the activation barriers for both clockwise and anti- 
clockwise directions are the same, as they are symmet- 
ric to each other. Activation barriers for translational 
concerted diffusion processes in all three directions (see 
Fig. 1) both for FF and HH dimers are reported in Ta- 
ble. I. Our results for concerted processes are 0.028 eV 
higher than the corresponding activation barriers for a 
dimer reported in Ref. 18. (This difference — as with 
those in what follows — may be due to the different inter- 
atomic potential employed in their study and ours.) 

Single-atom processes transform both FF and HH 
dimers into an FH-dimer. In this case one of the fee 



atoms in an FF-dimer or an hep atom in an HH-dimer 
moves to a nearest-neighbor hep or an fee site respec- 
tively, as shown in Fig. 2 (a)&(b) with the double-headed 
arrow. The activation barriers are 0.034 eV and 0.035 eV 
for hep and fee dimer, respectively. In the case of an FH- 
dimer, two types of single-atom diffusion processes are 
possible, as shown in Fig. 2: an fee atom moves to the 
nearest hep site in the direction of the open arrowhead, 
forming an HH-dimer, or an hep atom moves in the direc- 
tion of the solid arrowhead to the nearest fee site, forming 
an FF-dimer. The activation barriers for these processes 
are 0.028 eV and 0.024 eV, respectively. 



C. TVimer 

Depending on where a third atom is attached to the 
dimers shown in Fig. 2(a & b), there are four possible 
arrangements of atoms in a compact trimer: two types 
of fee timers - one centered around an hep site (F3H), 
the other centered around a top site (F3T) (see Figs. 3(a) 
& (d)) — and two types of hep trimers — one centered 
around an fee site (H3F), the other centered around a top 
site (H3T) (Figs. 3(c) & (b)). Although all four trimers 
have the same shape, their local environment is different, 
so that their adsorption energies are distinct, as are the 
activation barriers for their possible diffusion processes. 
F3T timer is the most energetically favorable: F3H, H3T 
and H3F are less energetically favorable by 0.006, 0.007 
and 0.0013 eV, respectively. It should also be noted that 
although trimers can take on non-compact shapes, the 
configurations depicted in Fig. 3 are the most frequently 
observed in our trimer simulations. 

In the case of F3T and H3T trimers, two types of 
concerted processes were observed, a non-diffusive con- 
certed rotation (clockwise and anti-clockwise) (Fig 3 (e) 
& (f)) and a diffusive concerted translation (in all three 
directions) (Fig 3 (d) & (b)). Concerted rotation pro- 
cesses transform an H3T timer into a F3T trimer and vice 
versa. The activation barrier for the concerted rotation 
processes for F3T trimer is 0.150 eV while that for those 
of the H3T trimer is 0.144 eV. Translation transforms an 
F3T timer into a H3F timer and vice versa. The acti- 
vation barrier for the concerted translations possible for 
these two trimers are 0.200 eV and 0.193 eV for F3T and 
H3T, respectively. The activation barrier of 0.200 eV for 
translational motion of F3T trimer is in agreement with 
the value reported in Ref 18. For F3H and H3F trimers, 
only concerted translation processes are possible; their 
activation barriers are 0.194 eV and 0.186 eV (the value 
reported for the same process in Ref 18 is 0.187 eV), 
respectively. Fig. 3 (a-b) & (c-d) reveal that these con- 
certed diffusion processes transform an F3H into an H3T 
timer and an H3F to an F3T trimer. Since the shape 
of these trimers is symmetric (see Fig 3 (a) & (c)), the 
activation barriers for their diffusion in all 3 possible di- 
rections are the same. 

As for single-atom processes in the case of a trimer: 
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FIG. 3. Possible arrangements of atoms in a trimer, with pos- 
sible concerted diffusion processes and their activation barri- 
ers, (a) — (d) Concerted translation: (a) F3H-all atoms on fee 
sites centered around an hep site; (b) H3T-all atoms on hep 
sites centered around a top site; (c) H3F-all atoms on hep 
sites centered around an fee site; (d) F3T-all atoms on fee 
sites centered around a top site, (e) & (f) Concerted rota- 
tion: F3T and H3T respectively. 



TABLE II. Activation barriers (in eV) for single-atom diffu- 
sion processes for an H3T compact trimer in the directions 
shown in Fig. 4. 



Type 


1 


2 


3 


4 


F3T 


0.439 


0.858 


0.858 


0.439 


F3H 


0.432 


0.875 


0.875 


0.432 


H3T 


0.436 


0.856 


0.856 


0.436 


H3F 


0.429 


0.872 


0.872 


0.429 



an atom can move in 4 different directions as shown in 
Fig. 4, resulting 2 different types of single-atom processes: 
directions 1 & 4 correspond to edge-diffusion processes 
which open up the trimer; directions 2 & 3 correspond to 
detachment processes (excluded from the present study, 
which is confined to diffusion of single whole islands, in 
which an island's integrity [and hence its size] is main- 
tained). We note that these processes move atoms from 
fee (hep) to nearest fee (hep) site. Activation barriers 
for processes in these 4 directions for different types of 
trimers are given in Table. II. Because these activation 
barriers are so high relative to those for concerted pro- 
cesses, single-atom processes were rarely observed in our 
simulations of trimer diffusion. 

We note that as island size increases, possible types 
of single-atom processes increases as well (though with 
the decamer, basically all possible types have appeared). 
Accordingly, it is convenient to defer detailed discussion 




FIG. 4. Single-atom processes possible for an H3T trimer. 
Activation barriers for the processes in these 4 directions for 
the 4 posible trimer configurations are given in Table. II 



of single-atom processes until later (Sub-section. IIIK) 



D. Tetramer 

Adding another atom to any of the trimers shown in 
Fig. 3 (a - d) results in the formation of a compact 
tetramer, diamond-shaped — with a long diagonal (along 
the line joining farthest atoms) and a short one (perpen- 
dicular to the long one), as shown in Fig. 5. Once again 
the fee island (Fig. 5(a)) is energetically more favorable 
than the hep one — in this case, by 0.009 eV. Three types 
of translational concerted diffusion processes are possible 
for each of the fee and hep tetramers, that is, one along 
each of the three directions specified in Fig.l. An exam- 
ple of a encerted fec-to-hep process (along direction 1) 
for a tetramer is shown in Fig. 5(a); its activation barrier 
is 0.213 eV. The reverse process (hep to fee) is shown 
in Fig. 5(b); its activation barrier is 0.204 eV (the value 
reported in Ref. 18 is 0.210 eV). Because the process 
in direction 3 is symmetric to that in direction 1, the 
energy barriers for these processes are identical, as are 
those for the reverse processes. The energy barrier along 
direction 2 is 0.313 eV from fee to hep and 0.304 eV from 
hep to fee. These values are systematically displayed in 
Table. III. 

Although the multi-atom processes shown in Fig. 5(c 
& d) have activation barriers lower than those of single- 
atom processes, they could not be found using the drag 
method. These processes thus had to be added man- 
ually to the database and their activation barriers ob- 
tained using the nudged elastic band method (NEB). In 
these multi-atom processes, two atoms move togehter in 
the same direction, the result is a shearing mechanism 
as shown in Fig. 5(c) & (d). For this shearing process, 
from fee to hep, the activation barrier is 0.285 eV; that 
for the reverse process from hep to fee is 0.276 eV. The 
drag method also finds single-atom processes, but be- 
cause in tetramers (as in islands of size 3 — 7) these have 
higher activation barriers than those of concerted pro- 
cesses, they were not observed during the simulations. 
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FIG. 5. Diffusion processes possible for a tetramer: (a) & 
(b) concerted diffusion along the short diagonal; (c) & (d) 
shearing processes. 



TABLE III. Activation barriers (eV) for the concerted 
tetramer translation processes shown in Fig. 5(a)&(b). 



Direction 


fee 


hep 


1 


0.213 


0.204 


2 


0.313 


0.304 


3 


0.213 


0.204 



E. Pentamer 




FIG. 6. Examples of concerted diffusion processes for a com- 
pact pentamer, along with their activation barriers. 

The compact shapes of a pentamer can be obtained by 
attaching an atom to a diamond-shaped tetramer. Al- 
though the geometries of compact pentamer clusters thus 
obtained are the same, the island's diffusion is crucially 
affected by where this additional atom is placed: attach- 



ment of an atom to an A-type step-edge of an fee tetramer 
results in the long A-type step-edge pentamer shown in 
Fig. 6(c); attachment of an atom to a B-type step-edge 
of the same tetramer results in the long B-type step-edge 
pentamer shown in Fig. 6(a)); the corresponding results 
of attaching an atom to an hep tetramer are shown in 
Figs. 6(b) and (d), respectively. The most energetically 
favorable of these is the fee pentamer with a long A-type 
step-edge (Fig. 6(c)); less favorlable by 0.005 eV are the 
fee pentamer with a long B-type step-edge (Fig. 6(a)), 
by 0.011 eV the hep pentamer with a long A-type step- 
edge, and by 0.017 eV the hep pentamer with a long 
B-type step-edge. That is: as usual, fee islands are more 
stable than hep ones. And, within each of those types, 
pentamers with a long A-type step-edge are more stable 
than those with a long B-type step-edge. 




FIG. 7. The single-atom processes that convert the long A- 
type step-edge pentamer to the long B-type step-edge pen- 
tamer. 



TABLE IV. Activation barriers (eV) without parantheses are 
for concerted-translation processes of pentamers with a long 
A-type step-edge, as shown in Figs. 6 (c)&(b); barriers in 
parentheses are for such processes for pentamers with a long 
B-type step-edge, as shown in Figs. 6 (a)&(d). 

Directions fee hep 

A (B) A (B) 

1 0.348 (0.301) 0.342 (0.284) 

2 0.348 (0.353) 0.342 (0.337) 

3 0.295 (0.353) 0.289 (0.337) 



In our simulations we found that compact pentamers 
diffuse mostly via concerted diffusion processes, which 
displace the island as a whole from fec-to-hep or vice- 
versa. Fig. 6 shows concerted diffusion processes along 
direction 1 for the long B-type step-edge pentamer and 
along direction 3 for the long A-type step-edge pentamer. 
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Table IV displays activation barriers for concerted pro- 
cesses in all 3 directions for both types of pentamer. 
Fig. 7 (a),(b) & (c) shows the single-atom processes that 
transform an fee pentamer from long A- type (= short B- 
type) to a short A-type (= long B-type) cluster, with the 
activation barrier for each. 



F. Hexamer 

Depending on whether a sixth atom is attached to 
a long A-type or to a long B-type step-edge pentamer, 
there are 3 possible compact shapes for a hexamer: (1) 
when an atom is added in such a way as to extend the 
shorter edge of either a long A-type or a long B-type 
step-edge pentamer, the result is one of the parallclo- 
gramic hexamers shown in Fig. 8; (2) when an atom is 
attached to the long edge of either type of pentamer, 
the result is the one of the irregular hexamers shown in 
Fig. 9; (3) when an atom is added to the shorter edge 
of either type of pentamer, the result is one of the trian- 
gular hexamers with all step edges of either the A-type 
(Figs. 10(b) &(c)) or of the B-type (Figs. 10(a) &(d)). 

Table. V shows the order of relative stabilties of the 




FIG. 8. Parallelogramic hexamers obtained by extending (by 
one atom) the shorter edge of either the long A-type or long 
B-type pentamers: (a) fee cluster; (b) hep cluster. The activa- 
tion barriers indicated are for concerted diffusion in direction 
1. 

hexamers most frequently obseved in our simulations. It 
reveals that, for a given shape, hexamers on fee sites are 
more energetically favored than those on hep sites, and 
that among hexamers on fee sites (as for those on hep 
sites), clusters in which A steps are longer than B steps 
are more stable. 

Figs. 8-10 also show concerted diffusion processes (in 
direction 1) for these hexamers, together with the acti- 
vation barriers for each. Tables. VI & VII give activation 
barriers for the hexamers shown in Figs. 8 & 9, respec- 
tively. Since triangular hexamers (Fig. 10) are symmet- 
ric, their activation barriers for concerted diffusion are 
same in all three directions. 

Fig. 11 shows the most frequently observed multi-atom 
processes for a hexamer - shearing processes in which a 
dimer moves along the A-type step-edge of the cluster 



TABLE V. Relative stabilities of the hexamers most fre- 
quently observed in our simulations. P = Parallelogram; I 



= Irregular; T 


= Triangular. 






Type 


Shape 


Description 


Energy (eV) 


Reference 


fee 


p 


equal A & B steps 





Fig: 8 fa) 


fee 


I 


edge atom on A step 


0.011 


Fig 9(b) 




I 


edge atom on B step 


0.011 


not shown 


hep 


P 


equal A & B steps 


0.014 


Fie 8(b) 


fee 


T 


all A steps 


0.019 


Fig. 10(c) 


hep 


I 


edge atom on B step 


0.024 


Fig. 9(b) 


hep 


I 


edge atom on A step 


0.024 


not shown 


fee 


T 


all B steps 


0.029 


Fig. 10(a) 


hep 


T 


all A steps 


0.032 


Fig. 10(b) 


hep 


T 


all B steps 


0.042 


Fig. 10(d) 




FIG. 9. Irregular hexamers obtained by attaching an atom to 
the long edge of a pentamer: (a) fee cluster; (b) hep cluster. 
The activation barriers indicated are for concerted diffusion 
in direction 1. 



from sites of one type to the nearest-neighbor sites of the 
same type. Figs. 11 (a) & (b) show this kind of diffusion 
process for an hep cluster and Figs. 11 (c) & (d) for an fee 
cluster. Although this dimer shearing process does not 
much displace the center of mass of a hexamer; it does 
have a striking consequence: it converts a parallelogramic 
hexamer (Fig. 8) into an irregular hexamer, (Fig. 9) and 
vice- versa. 



G. Heptamer 

On an fcc(lll) surface, an heptamer has a compact 
closed-shell structure with each edge atom having at least 
three nearest-neighbor bonds, as shown in Fig. 12. Our 
SLKMC simulations (keep in mind here the range of 
temperatures to which they were confined) found that 
heptamer diffuses exclusively via concerted diffusion pro- 
cesses, which displace the cluster from fec-to-hep and vice 
versa; the barriers for which are shown in Fig. 12. That 
these processes will predominate can also be concluded 
from the fact that the effective energy barrier for hep- 
tamer diffusion (cf. Table. 23) is close to the average of 
the activation barriers shown in Fig. 12. 

Since the compact heptamer has a symmetric shape, 
activation barriers in all three directions are the same as 




FIG. 10. Triangular hexamers obtained by adding an atom 
to the short edge of a pentamer: (a) fee hexamer with B- 
type step edges; (b) hep hexamer with A-type step edges; (c) 
fee hexamer with A-type step edges; (d) hep hexamer with B- 
type step edges. The activation barriers here are for concerted 
diffusion in direction 1. 

TABLE VI. Activation barriers (eV) of the concerted trans- 
lation processes in all three directions for the hexamer shown 
in Fig.8 



Directions 


fee 


hep 


1 


0.374 


0.360 


2 


0.466 


0.451 


3 


0.254 


0.240 



those shown in Fig. 12. Again, the fee island is more 
energetically favorable than its hep counterpart - in this 
case by 0.015 eV. 

H. Octamer 

Compact octamers have two distinct orientations, one 
with two long A-type step-edges, the other with two long 
B-type step-edges, as shown in Fig. 13 (a) & (b) for an fee 
octamer. Octamers with long A-type step-edges (Fig. 13 
(a)) can be obtained by attaching an atom to any Pi- 
type step-edge of a compact heptamer, while a compact 
octamer with long B-type step-edges results from attach- 
ing one to any A-type step-edge. Again: the fee islands 
are more energetically favorable than the hep ones, and 
within each type, islands with long A-type step edges are 
more stable than those with long B-type step-edges. 

A compact octamer diffuses via concerted diffusion 
processes, as shown in Fig 14. The activation barrier of 
a concerted diffusion process depends on whether the oc- 



TABLE VII. Activation barriers (eV) of concerted transla- 
tions processes in all 3 directions for the hexamers shown in 
Fig.9 



Directions 


fee 


hep 


1 


0.393 


0.380 


2 


0.397 


0.383 


3 


0.391 


0.378 




FIG. 11. Dimer shearing processes in case of a hexamer along 
with their activation barriers: (a) & (b) all-hep hexamers; (c) 
& (d) all-fee clusters. 



tamer has long A- or long B-type step-edges. As Fig. 14 
shows, a concerted diffusion process converts a long A- 
type step-edge fee octamer into a long B-type step-edge 
hep cluster, and vice-versa. Table. VIII reports the acti- 
vation barriers for concerted diffusion processes in all 3 
directions for both orientations for an fee as well as an 
hep (see Fig. 13) octamer. Concerted diffusion processes 
in directions 2 and 3 are the most frequently observed 
processes in octamer diffusion. 

Although an octamer diffuses primarily via concerted 
processes we found in our simulations that both multi- 




FIG. 12. Concerted diffusion processes and their activation 
barriers in direction 1 for a heptamer. 
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FIG. 13. Possible orientations for a compact fee octamer: (a) 
with long A-type step edges; (b) with long B-type step edges. 




FIG. 14. Examples of concerted diffusion processes in di- 
rection 2 for an octamer: (a) compact fee octamer with long 
A-type step-edges; (b) compact hep octamer with long B-type 
step-edges. 



atom and single-atom processes are also relatively com- 
mon. As mentioned before, we defer comprehensive dis- 
cussion of single-atom processes to section IIIK. Here 
we discuss multi-atom processes particular to octamers. 

Multi-atom processes involving shearing and reptation 
are shown in Figs. 15 & 16, respectively. In shearing 
processes, part of the island (more than one atom) moves 
from fee to the nearest fee sites, if all the island is initially 
on fee sites - or from hep to hep sites, if all of it initially 
sits on hep sites. Figs. 15 (a) & (b) show trimer shearing 
processes within an octamer, along with their activation 
barriers, while Figs. 15 (c) & (d) show dimer shearing 
processes, with their activation barriers. 

Reptation is a 2-step diffusion process. In the case 
of an fee island, the entire island diffuses from fee to 
nearest-neighbor hep sites in two steps. In the first, part 
of the island moves from fee to nearest-neighbor hep sites, 
leaving part of the island on fee sites and part on hep 



TABLE VIII. Activation barriers (in eV) of the twelve con- 
certed diffusion processes for compact octamers. 

Directions fee hep 

A (B) A(B) 

1 0.589 (0.585) 0.567 (0.571) 

2 0.491 (0.484) 0.468 (0.469) 

3 0.491 (0.484) 0.468 (0.469) 




FIG. 15. Example of shearing diffusion processes within an 
octamer along with their activation barriers. 



(a) (b) 




FIG. 16. Successive sub-processes (or steps) involved in an 
octamer reptation diffusion mechanism. 



sites. In the next, the remainder of the island initially 
on fee sites moves to hep sites. Figs. 16 (a)-(d) shows 
various steps (subprocesses) of a reptation process, with 
their activation barriers. 



I. Nonamer 

For a nonamer, we observed all types of diffusion pro- 
cesses - single-atom, multi-atom and concerted. The 
most frequently observed were two types of single-atom 
mechanisms: edge-diffusion processes along an A- or a B- 
type step-edge and corner rounding (Fig. 21). The non- 
amer is the smallest island for which concerted processes 
are not the most frequently picked (the concerted pro- 
cesses for the most frequently observed nonamer configu- 
rations [compact or nearly so] are shown in Fig. 17). Even 
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FIG. 17. Concerted diffusion processes and their activation 
barriers for nonamers. 




FIG. 18. Dimer shearing processes and their activation bar- 
riers for compact nonamers. 



so, concerted processes contribute the most to island dif- 
fusion: that is, the displacement they produce in the non- 
amer's center of mass is far greater than that produced by 
single-atom processes, despite the far greater frequency 
of the latter. This is reflected in the fact that the effective 
energy barrier for nonamer (cf. Table XII) is much closer 
to the average activation barrier for concerted processes 
(cf. Table. IX) than for that of single-atom processes 
(cf. Table IX) . The fact that the effective activation bar- 
rier is slightly higher than the average energy barrier for 



concerted processes is due mainly to the contribution of 
kink processes, which do contribute somewhat to island 
diffusion. 

TABLE IX. Activation barriers (eV) of concerted transla- 
tions processes in all 3 directions for the nonamers shown 
in Figs.l7(a-d). 



Directions 


fee 


hep 


fee 


hep 




(a) 


(b) 


(c) 


(d) 


1 


0.520 


0.499 


0.486 


0.465 


2 


0.626 


0.605 


0.486 


0.465 


3 


0.605 


0.583 


0.693 


0.672 



The most frequently observed multi-atom processes 
are the four forms of dimer shearing along an A-type 
step-edge shown in Fig. 18(c)-(f), as have been discussed 
above for island of sizes 6 & 8. The activation barriers 
for these dimer shearing processes are lower than those 
for single-atom diffusion processes along an edge and also 
for some corner rounding processes. Reptation processes 
also show up, but only when the nonamer is non-compact 
(we do not illustrate these here) 37 . 

J. Decamer 

Even in the case of a decamer, we have observed 
single-atom, multi-atom and concerted diffusion pro- 
cesses. Single-atom diffusion processes are the most fre- 
quently observed. The most frequently observed compact 
shape of decamer during our simulations is that shape 
shown in Fig. 19, which has the same number of A- and 
B-type step edges. As usual, an fee cluster is energeti- 
cally more favorable than an hep cluster. 

TABLE X. Activation barriers (eV) of concerted diffusion pro- 
cesses in all 3 directions of decamer as shown in Fig. 19 (a) & 
(b). 



Directions 


fee 


hep 


1 


0.661 


0.638 


2 


0.700 


0.677 


3 


0.700 


0.677 



For the shape shown in Fig. 19, the most frequently 
observed concerted diffusion processes are those shown 
in the same figure, along with their activation barriers 
reported in Table X. It can be seen from the Table 23, 
that the effective energy for decamer diffusion is close 
to that of the average energy barrier of these concerted 
processes. That is why decamer diffusion is dominated 
by concerted processes. 

As with the nonamer, a decamer also undergoes multi- 
atom processes (shearing and reptation). Of these, the 
most frequently observed is dimer shearing along an A- 
type step-edge, similar to what has been discussed for 
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FIG. 19. Frequent concerted diffusion processes and their 
activation barriers for compact decamers. 




FIG. 20. Diffusion steps (sub-processes) in decamer reptation. 
Note that the decamers in (b) and (c) are identical, though 
the arrows indicate different processes. 

clusters of size 6, 8 and 9. Fig. 20 shows the sub-processes 
in a reptation process, along with the activation barrier 
of each. 



K. Single-atom Processes 

In this section we provide detail about single-atom pro- 
cesses: edge-diffusion, corner rounding, kink attachment, 
and kink detachment, as shown in Fig. 21 for an hep is- 
land. Their corresponding activation barriers and those 
for their fee analogues are given in Table XL In each 
single-atom process, an atom on an fee site moves to a 
nearest-neighbor fee site, while an atom on hep site moves 
to a nearest-neighbor hep site. The activation barriers for 
single-atom processes depend not only on whether the 
atom is part of an fee sland or an hep island, but also on 
whether the diffusing atom is on an A-type or a B-type 



2 2 4 4 




FIG. 21. Single-atom processes for an hep island (Though 
analogous processes occur for an fee island, we do not illus- 
trate them here). The index numbers designate the processes 
described in Table XI, which gives the activation barriers for 
each. 



step-edge. 

In classifying single-atom processes in Table XI we 
have used the notation X ni U — > Y„,V, where where X 
or Y = A (for an A-type step-edge) or B (for a B-type 
step-edge) or K (for kink) or C (for corner) or M (for 
monomer); n^ = the number of nearest-neighbors of the 
diffusing atom before the process; n/ = the number of 
that atom's nearest neighbors after the process. U or V 
= A or B (for corner or kink processes) or null (for all 
other other process types). 

For example, process 1, B 2 — > B 2 , is a single-atom 
B-step edge process in which the diffusing atom has 2 
nearest-neighbors before and after the process. Process 
3, CiB —¥ B2, is a corner rounding process towards a 
B-step, the diffusing atom starting on the corner of a 
B-step with one nearest-neighbor and ending up on the 
B-step with two nearest-neighbors. In process 10, C 2 A 
— > CiB, the diffusing atom begins on the corner of an 
A-step having two nearest-neighbors and ends up on the 
corner of a B-step with only one nearest-neighbor. 
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TABLE XI. Activation barriers (eV) of single-atom processes 
for both fee and hep islands. The index numbers refer to 
the types of processes illustrated in Fig. 21. See text for 
explaination of the notation used to classify the process types. 

Index no. Process type fee hep 



1 


B 2 -> 


B 2 


0.454 


0.448 


2 


B 2 -» 


M 


0.821 


0.815 


3 


CiB 


-> B 2 


0.177 


0.173 


1 


CiB 


-> M 


0.458 


0.455 


5 


CiB 


-> A 2 


0.040 


0.038 


6 


C]B 


-> CiB 


0.540 


0.585 


7 


CiB 


-> M 


0.811 


0.809 


8 


A 2 -» 


M 


0.795 


0.794 


9 


A 2 -» 


A 2 


0.326 


0.307 


10 


C 2 A 


-> CiB 


0.399 


0.397 


11 


C 2 A 


-» M 


0.787 


0.785 


12 


B 2 -> 


K 3 B 


0.415 


0.298 


13 


K 3 A 


-> A 2 


0.601 


0.701 


14 


A 2 -» 


K 3 A 


0.302 


0.389 


15 


K 3 B 


-> B 2 


0.729 


0.597 


16 


K 3 B 


-» CiB 


0.731 


0.787 


17 


K 3 B 


-> M 


1.138 


1.150 


18 


K 3 B 


-> CiB 


0.820 


0.759 



IV. DIFFUSION COEFFICIENTS AND 
EFFECTIVE ENERGY BARRIERS 

We start our SLKMC simulations with an empty 
database. Every time a new configuration (or neighbor- 
hood) is turned up, SLKMC-II finds on the fly all possible 
processes using the drag method, calculates their acti- 
vation barriers and stores them in the database as the 
simulation proceeds. Calculation of energetics occurs at 
each KMC step during initial stages of the simulation 
when the database is empty or nearly so, and ever less 
frequently later on. Recall that the types of processes 
and their activation barriers are dependent on island size, 
each one of which requires a separate database that can- 
not be derived from that for islands of other sizes. 

We carried out 10 7 KMC steps for each island size at 
temperatures 300K, 400K, 500K, 600K and 700K. We 
calculated the diffusion coefficient of an island of a given 
size using Einstein Equation: 44 D = lim f _j.oo (Rcm {t) — 
-Rcm(O)] 2 ) /2dt, where Rcniit) is the position of the cen- 
ter of mass of the island at time t, and d is the dimen- 
sionality of the system, which in our case is 2. Diffusion 
coefficients thus obtained for island sizes 1 — 10 at var- 
ious temperatures are summarized in Table XII. At 300 
K, diffusion coefficients rang e from 1.63 x 10 11 A 2 /s for 
a monomer to 8.66 x 10 01 A 2 /s for a decamer. Effective 
energy barriers for islands are extracted from their re- 
spective Arrhenius plots (Fig. 22) and also summarized 
in Table XII. Fig. 23 plots effective energy barrier as a 
function of island size. It can be seen that the effective 




1/k T [eV]" 1 

B 

FIG. 22. Arrhenious plots for 1—10 atom islands, 
energy barrier increases almost linearly with island size. 




Island size (atoms) 



FIG. 23. Effective energy barriers of 1 — 10 atom islands as a 
function of island size. 



V. FURTHER DISCUSSION AND 
CONCLUSIONS 

To summarize: we have performed a systematic study 
of the diffusion of small Ni islands (1-10 atoms) on 
Ni(lll), using a self-learning KMC method with a newly- 
developed pattern recognition scheme (SLKMC-II) in 
which the system is allowed to evolve through mech- 
anisms of its choice on the basis of a self-generated 
database of single-atom, multiple-atom and concerted 
diffusion processes (each with its particular activation 
barrier) involving fcc-fcc, fec-hep and hep-hep jumps. 
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TABLE XII. Diffusion coefficients ((A 2 /s) at various temperatures and effective energy barriers for Ni islands. 



Island Size 



300K 



400K 



500K 



600K 



700K E eff (eV) 



1 


1.63 


X 


10 11 


2.85 x 10 11 


3.99 x 10 11 


5.00 


X 


10 11 


5.87 


X 


10 11 


0.058 


2 


8.14 


X 


10 09 


2.88 


X 


10 10 


6.54 


X 


10 10 


1.13 


X 


10 11 


1.64 


X 


10 11 


0.136 


3 


1.09 


X 


10 09 


7.42 


X 


10 09 


2.30 


X 


10 10 


5.00 


X 


10 10 


8.23 


X 


10 10 


0.196 


4 


7.64 


X 


10 08 


5.83 


X 


10 09 


2.02 


X 


10 10 


4.91 


X 


10 10 


9.93 


X 


10 10 


0.217 


5 


2.90 


X 


10 06 


8.42 


X 


10° 7 


6.59 


X 


10 08 


2.55 


X 


10 09 


7.36 


X 


10 09 


0.353 


6 


1.08 


X 


10 06 


7.39 


X 


10 07 


5.66 


X 


10 08 


2.69 


X 


10 09 


8.12 


X 


1Q 09 


0.400 


7 


6.24 


X 


10 04 


6.40 


X 


10 06 


1.01 


X 


10 08 


6.37 


X 


10 08 


2.37 


X 


10 09 


0.477 


8 


1.38 


X 


10 04 


1.50 


X 


10 06 


2.31 


X 


10 07 


1.57 


X 


10 08 


5.91 


X 


10 08 


0.482 


9 


5.48 


X 


1Q 02 


8.77 


X 


1Q 04 


1.89 


X 


10 06 


1.35 


X 


10 07 


6.23 


X 


10 07 


0.525 


10 


8.66 


X 


10 01 


2.25 


X 


10 04 


7.20 


X 


10 05 


7.97 


X 


10 oe 


4.95 


X 


10 07 


0.627 



We find that concerted diffusion processes contribute the 
most to the displacement of the center of mass (i.e., to 
island diffusion), while single-atom processes contribute 
the least. As for multi-atom processes (reptation or 
shearing): while these produce more displacement than 
the latter, they are hardly ever selected. Though the 
energy barriers for reptation processes are small com- 
pared to those for concerted diffusion processes, reptation 
occurs only when an island is transformed into a non- 
compact shape, as happens only rarely in the temprature 
range to which our study is confined. In contrast, though 
shearing occurs with close-to-compact shapes (which ap- 
pear more frequently than do non-compact shapes, but 
with islands of certain sizes), the barriers for these pro- 
cesses are higher than those for reptation. Finally, al- 
though for all island sizes, island diffusion is primar- 
ily dominated by concerted diffusion processes, the fre- 
quency of occurrence of both single-atom and multi-atom 
processes does increase with increase in island size, owing 
to increase in the activation barrier for concerted diffu- 
sion processes with island size. 

Fig. 24 shows the number of each type of diffusion pro- 
cess (single-, multi-atom and concerted) collected during 
our SLKMC-II simulations for each island size (1-10), to- 
gether with the total number of processes of all types for 
each island size. (For the sake of clarity, the insert shows 
a log-linear plot of these quantities for island sizes 1-6). 
It can be seen that the number of processes accumulated 
increases with island size, and significantly so beyond 
the tetramer. It can also be seen from Fig. 24 that the 
overall increase in number of processes with island size 
is constituted predominantly by significant increases in 
single-atom processes and (to a lesser degree) multi-atom 
processes. Meanwhile the number of concerted processes 
accumulated in the database increases at a much slower 
pace with island size. 

This significant increase in single-atom processes is 
mainly due to the use of 10 rings to identify the neigh- 
borhood of an atom. Use of 10 rings corresponds to in- 
clusion of 5 nearest-neighbor fcc-fcc or hcp-hcp interac- 
tions. Elsewhere we show that 6 rings (which corresponds 
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FIG. 24. Distribution of single-atom, multi-atom, concerted 
and total processes for 1-10 atom islands accumulated in the 
database during SLKMC simulations. Inset shows the log- 
linear plot for up to the 6-atom island. 



to 3 nearest-neighbor interactions) offer a range of inter- 
action sufficient for accurately calculating the activation 
barriers for single-atom processes. 45 But it is essential to 
include the long-range interaction (and hence 10 rings) 
if one aims to accurately take into account multi-atom 
and concerted processes, the latter of which predominate 
in small-island diffusion. This significant increase in the 
number of processes with island size also justifies resort- 
ing to an automatic way of finding all the possible pro- 
cesses during simulations instead of using a fixed (and 
thus necessarily preconceived) list of events. 

As mentioned earlier, with increasing island size not 
only does the number of accumulated single-atom pro- 
cesses and multi-atom processes increase but also their 
frequency of occurrence. Still island diffusion is primar- 
ily due to concerted diffusion processes, since it is these 
that produce largest displacement of center of mass. This 
can be easily observed by comparing, for each island size, 
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the effective diffusion barriers given in Table XII with 
the activation barriers in the tables given in Sect. Ill for 
concerted diffusion processes: effective diffusion barriers 
more or less closely follow activation barriers for con- 
certed diffusion processes - except for the 9-atom island, 
in which the contribution of single-atom processes to the 
island's diffusion is significantly larger than for other is- 
land sizes. 

In conclusion, Ni small-island diffusion on Ni (111) is 
primarily due to to concerted diffusion processes, even 
though the frequency of their occurrence decreases with 



increase in island size. 
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